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Abstract 

The universe is viewed as a dust gas filling a sphere and floating in infinite empty space. Einstein's 
gravitational equations are applied to this case together with appropriate boundary values. The 
equations are solved for initial conditions chosen so as to describe the observed Hubble diagram. We 
find that the solution is not unique so that more astronomical observations are needed. However, those 
solutions which were found do not exhibit an accelerated expansion of the universe, nor - obviously 
then - do they need the notion of a dark energy driving such an expansion. 

We present this study as an alternative to the prevailing Robertson-Walker cosmology. 


1. Introduction 

The FRW model 

The currently popular cosmologjQ is much influenced by the need to understand or describe the 
observed relation between the apparent magnitude m of Type la supernovae and the redshift z of 
the light received from them. The observations, reported by Riess et al., 2, Amanullah et al. [ 4 J, 
Hicken et al. [§] and others, compiled in the Union 2.0 and 2.1 catalogues 6], are reproduced in the 
diagram of Fig. [I] as dots. The diagram is often called the Hubble diagram, because in some rough 
manner m determines the distance of the light source and z represents its speed; and Hubble Q was 
first to conjecture a relation between speed and distance. That conjecture gave rise to the notion 
of an expanding universfl The need to understand the (m,z)-correlation has led cosmologists to a 
revival of the cosmological constant^ A, - interpreted as a hypothetical repulsive "dark energy" -, 
and the prediction of an accelerated expansion of the universe. These notions are embedded in the 
concept of a four-dimensional Euclidean space with a three-dimensional surface, either spherical or 
flat or hyperspherical, and endowed with a Robertson-Walker metric. The three-dimensional surface is 
supposed to be homogeneously filled with matter, a strongly simplifying assumption that reduces the 
Einstein equations of gravitation to a simple set of ordinary differential equations. Friedman 0,112 
has solved the equations and accordingly this theory is often called the FRW-cosmology. The study of 
the cosmic radiation background has led FRW-cosmologists to the assumption that the universe is in 
fact flat. And if this is the case, the Hubble diagram requires that the dark energy makes up ca. 70% 
of the mass of the universe (!) (see [l|, p.48). 

Present model 

In the present paper we investigate a less arcane model: We view the universe as a sphere of 
the time- dependent radius R(t) filled with matter and floating freely in infinite empty spac43- The 
matter is supposed to be distributed isotropically with respect to the center of the sphere and it is 


* Corresponding author 

Email addresses: ingo.mueller@alumni.tu-berlin.de (Ingo Muller), wolf.weiss@alumni.tu-berlin.de (Wolf 
Weiss) 
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Figure 1: Dots: Observations of distance moduli p = m — M and redshifts 2 for Type la supernovae. M is the absolute 
magnitude; it is the same for all Type la supernovae and it has the value -19. 

Graph: The (//. 2)- plot for our model of the universe fits well to observed values. 


considered as a dust gas so that pressure and internal energy are negligible and temperature plays no 
role. Cosmic radiation is neglected as well and thus the only agent dictating the motion of the dust 
gas is the gravitational effect on it. 

The motion is governed by the Einstein equations of gravitation, - without the cosmological con¬ 
stant A. We solve these equations for appropriate boundary conditions, - in the center, at infinity, 
and at R(t) -, and for appropriate initial conditions at some time ti in the past, when the expansion 
observed by Hubble is already in progress. The initial conditions are choosen so as to provide a good 
fit between the calculated /i(z)-curve and the observed data, see Fig. [TJ the method is known in math¬ 
ematics as the solution of an inverse problem. From the initial time onwards we follow the continued 
expansion for t > ti and, calculating backwards for t < ti, we see that the universe has emerged from 
a sphere with the Schwarzschild radius. 

We come to the conclusion that the observed p,(z )~data are compatible with different solutions of 
the Einstein equations,- always without A. That is to say that different initial data lead to the same 
/i(z)-curve. We exhibit two solutions explicitly and none of them shows an accelerated expansion nor 
do they require dark energy. The non-uniqeness of the solution of our inverse problem means that more 
astronomical observations are required to find a unique solution. Which additional observations are 
feasible? And which solutions of the Einstein equations describe them well? Those questions represent 
the cosmological conundrum. 

In our paper we choose Schwarzschild coordinates ( t , r, i9, ip) for calculations , i.e. for the formulation 
and reformulation and solution of the Einstein equations. In these coordinates we have a fairly simple 
diagonal metric tensor with only the radial and the temporal component as initially unknown functions 
of t and r, while the angular components are like those in Euclidean space. For the interpretation of 
the results we replace the Schwarzschild coordinate time t by the proper time of the observer in the 
center of the sphere. 

It is true that in our model of the universe we, the observers of luminosities and redshifts, are 
situated in the center of the spherical dust gas or close to it. This violates the cosmological principle , 
according to which we should not occupy a privileged place in the cosmos. We do accept that violation 
because we believe that the privileged place is less counter-intuitive than is a hypothetical dark energy. 
Besides in Chapter 6 we conjecture that the position of our observer need not necessarily be restricted 
to the central point of the sphere, rather it may he within a fairly extensive sphere. 

2. Einstein equations 

2.1. General form 

The Einstein equations read (without cosmological constant) 

Gap = 8 n^T a0 , ( 1 ) 
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where G a p = R a p — —Rg a g is the Einstein tensor and T a p is the energy-momentum tensor of the 

matter. 7 is the gravitational constant. R a g is the Ricci tensor and R = g a ^ R a g is the Ricci scalar. 
The Ricci tensor is related to the Christoffel symbols T^ which in turn are related to the metric tensor 
9y, v - We havc0 
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2.2. Metric tensor. Speed and trajectory of light. 

For the present purpose, - the consideration of a radially expanding sphere, isotropically filled with 
matter - we assume that the metric tensor has the form 


9ij = 


S 


r sm' 


0 J 


(3) 


Space-time coordinates (t, r , $, <p) for which the metric tensor has this form are called Schwarzschild 
coordinates. Oppenheimer & Volkoff 141 and Oppenheimer & Snyder [15] have used Schwarzschild 
coordinates for the description of a neutron star. Here we use them for the description of the universe 
as a whole: A spherical dust gas in infinite empty space. A sphere of radius r in Schwarzschild 
coordinates obviously has the surface area 4nr 2 just like in Euclidean space, see also Misner et al [13!]. 
The components Z and S may be functions of t and r, and we rely on the Einstein equations to find 
the dependence. Z(t,r ) must be negative, because without gravitation Z must be equal to -1. S must 
be equal to 1 in that case. 

Since Z(t,r) and S(t , 7 ') are components of the metric tensor, we feel justified to assume that they 
must be continuous and smooth functions, apart possibly from singular points. 

It follows that the infinitesimal distance element ds 2 = g a gdx a dx ^ in space-time is given in terms 
of the coordinate increments dt and dr 


ds 2 = Zc 2 dt 2 + Sdr 2 . (4) 

In a local and momentary Lorentz frame at rest with a particle of matter we denote the coordinate 
increments of time and radial distance by dO and dg so that we have 

ds 2 = — c 2 dd 2 + dg 2 . (5) 

dO and dg are called the increments of proper time and proper distance associated with the particle. 
Hence follows a relation between dt and dO for fixed g 


_ / T/2 c (g r 

d6 = VMl/l - whereC = [- 


( 6 ) 


is the velocity of a material particle in the (t, 7 ')-system. Similarly a relation between dr and the 
increment dg of proper distance at fixed 8 readfl 


de= VsJl-^dr. 


(7) 


5 Greek indices run from 0 to 3, Latin ones from 1 to 3, such that x° = ct, x b = (r, <£>). is the metric tensor in 

space-time. g a p and its inverse g a @ may be used to lower and raise indices in the usual manner. 

6 We drop the angular part of ds 2 for brevity; it is unimportant for radial motion. 

7 While j6j is an obvious consequence of ((4j), (j5j, the relation J7]) requires a few intermediate lines of calculation with 

f dt\ VS 

partial derivatives in order to determine that — = —- holds. 
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The trajectory of a radial light ray in the (f, r)-system is governed by the equation ds 2 = 0 so that, 
by (U), we have for the velocity c x of light moving radially 



\c x \ denotes the speed of light in the (t, r)-system. Obviously, from ©, Q) the speed V of a material 
particle is bound by the speed of light in the (t, r)-system 


\V\ < c x 



(9) 


Later we shall be interested in the trajectory of a light ray that reaches the observer of our spherical 
universe in the center of the sphere at the present time t = 0, i.e. at the event (t,r) = (0,0). That 
trajectory - denoted by rx(t) - must be obtained by integration from (|5j). Naturally we need explicit 
functions Z(t,r) and S(t,r) for the purpose. Therefore such a calculation has to be postponed until 
we have obtained those functions from the Einstein equations. 


2.3. Specific form of the Einstein tensor 

It is a cumbersome task to calculate the Christoffel symbols, hence the Ricci tensor, and hence the 
Einstein tensor from the specific form J3|) of g a0 . However, when the calculations are done, it turns 
out that the Einstein tensor contains only four essential non-zero components and they read 
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()' and ()* denote partial derivatives with respect to r and t respectively. G vv is not zero, but it is 
equal to G&& to within a factor r 2 sin 2 d and may therefore be ignored, because the components 
and Tptp of the energy-momentum tensor both vanish in our case as we shall presently see. 


2-4- Energy-momentum tensor for a dust gas 

We assume that the matter in the sphere is a dust gas. This means that there are no viscosities 
and, in addition, that the pressure may be neglected. If that is true, the energy-momentum tensor is 
given by 

T a0 = aU a U 0 , (11) 


where a is a scalar rest-mass density of the matter. The idea is, of course, that the matter in the 
universe is "smeared out" to form a continuum; each point is then occupied by a "particle" of matter 
in the sense of continuum mechanics. 

U a are the covariant components of the four-velocity and U a = g a/3 U /3 are its contravariant com¬ 
ponents. For radial expansion we have 
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V = f — I is the radial velocity of the gas as before, and 9 is the proper-time. We use © and © 
to write (fl2l) in the form 
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dr 

Note that the r-component ( — J of the 4-velocity equals 
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V 


v'M./TTw 


( dv \ 

so that | — | varies between 0 and oo when \V\ varies between 0 and \c x \, the speed of light. 

V Q J Q 

By dill) , d!3D the energy-momentum tensor assumes the form 
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2.5. Specific form of the Einstein equations 

Elimination of G a p and T a p between ©, m, and (fill) gives the specific form of the Einstein 
equations for the radial expansion of an isotropic sphere filled by a dust gas. We obtain 
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We rely on these four equations for the determination of the four fields Z(t,r), S(t,r), a(t,r), and 
V (f, r). First, however, we rewrite these equations aiming for a set of equations that can be recognized 
as a relativistic generalization of the Newtonian non relativistic - set appropriate for a self-gravitating 
dust gas. 

First of all we note that mi may be written in the form 


r ,!- 1 


= -2 V 


or, by integration from 0 to r with m(t,r) = / m'(t,a)dc 
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m is the mass within a sphere of radius r. Indeed, we observe that the integrability condition for S 
implied by (fl5l) i 9 reads 

m' + (m'V)' = 0. (17) 

It follows that m! is the density of a conserved quantity. Therefore m! is the mass density of a spherical 
shell referred to the thickness dr of the shell. \m! must not be confused with the scalar rest-mass density 
which is denoted by a, cf. hum Thus, by (fThl) the dependent field S(t,r ) may be replaced by the 
more intuitively appealing field m(t,r). This field satisfies a very simple differential equation which 
we derive as follows. We observe that GHDi multiplied by V and added to (fThlh provides the relation 


In ■ 


S 


S- 1 


V 


In ■ 


S- 1 


V 

= -, or by GHD 


m 

In — 


V 


m 

In — 


=-, or, finally 
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m + Vm' = 0 . 


It follows that the mass inside a spherical surface of radius r is constant for the observer moving with 
that surface. That result is eminently plausible, of course, and we write it as 


dm n , d() ... Tr/N/ 
— = 0 , where — = () + VQ . 


(18) 


The formidable equation (11511 a may be given a somewhat more intuitively appealing form by elim¬ 
ination of S, Z", S , S' and by use of (HU), CGD. We obtain after some calculation 
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This equation may be seen as the relativistic generalization of the equation of motion of a self-gravitating dust gas. 
Indeed, in the non-relativistic case, i.e. when Z = — 1 and V 2 /c 2 « 0, we obtain from m 

dV m 
~di ~ 

which is the r-component of the non-relativistic equation of motion with the classical attractive gravitational acceler¬ 
ation on the right-hand-side. If an expansion is in progress, so that V > 0 for all r > 0 holds, that term decelerates 
the expansion. Note however, that in oi the second term on the right-hand-side is positive, i.e. that term causes an 
accelerated expansion. In Section 5.5 the effect of that term is made explicit after the Einstein equations have been 
solved and the fields Z(t,r), m(t,r), cr(t, r), and V(t,r) have been determined. Note also that the second term on the 
right-hand-side of 01 becomes singular when the sphere has the radius rg = 2'yM/c 2 which is the Schwarzschild radius 
for the total mass M. 


Finally we may eliminate the acceleration dV/dt from (flTH) by use of (flSlh . 
Making use of (TTTjl . (fl8l) we thus obtain 
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(ar 2 )'+ (ar 2 V) = 

c 2 r 


( 21 ) 


This is clearly the relativistic generalization of the non-relativistic law of mass-conservation, since the right-hand-side 
is negligible in that case. Note that non-relativistically there is no essential difference between ar 2 and the mass density 
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In summary we may now say that so far the original four Einstein equations are replaced by the 
four equations 
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Alternatively - replacing 


1 by its corollary (l20l) - we may write this set of equations as 
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Yet another arrangement of these equations is as follows 
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This results from replacing ( 1221) 4 by its corollary KUl) . The three sets (1221) through (1M1) are equivalent; 
their different forms are needed here for arguments of interpretation and in the numerical solution of 
the set. 


2.6. Non-dimensional variables 

For the subsequent numerical calculations we introduce dimensionless dependent and independent 
variables based on the speed of light c, the total mass M of the universe and its present radius I?(0). 
Note that both M and -R(O) are a priori unknown; they must be chosen so as to be compatible with 
the p.(z)- curve. We define 


r, t = 


ct 
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_ 47T 
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Also we introduce the dimensionless parameter 
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which represents the Schwarzschild radius of the universe referred to its present radius f?(0). 


(25) 


(26) 


3. External solution for S(t,r ) and Z(t,r) 

In the space outside the matter-filled sphere the rest mass density a is equal to zero so that mi 
may be integrated to give S = K{t)/{ 1 — K(t)/r). K(t) is a constant of integration whose value is 
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given by © as (2 7 /c 2 ) M , because at the surface of the sphere S must be continuous and equal to 
the value given by (fTTH) with m = M = m(t, R(t)) which is the constant total mass. Thus we have for 
the external solution of S 


SfiT.t - 


1 27 M 

J- „ 


forr > R(t), ot Sext = 


I-Q7 


for f > R{t). 


(27) 


If we require S(t, r ) to be continuous and smooth at the surface, we must obviously have m'(t, R(t)) = 0. 
Z ex t, the outer solution of Z(t,r), follows by insertion of (12711 into ( 1221 ) 9 . and integration: Z ext = 
( 2y M \ 

k(t) 1 - 5 -, where kit) is a constant of integration whose value must be equal to — 1 , so that 

V c r ) 

Z —> — 1 for r —> 00 , where gravitation is negligible. Thus we have for the outer solution of Z 

21 - . 

Zext = -1 + —— forr > R(t) or Z ext = -1 + Q- forr > R(t). (28) 

c r r 

Neither S ex t nor Z ex t depend on time. And obviously S ex t is singular at 

2y 

r =r s = -4-M, or r = r s = Q, (29) 

c z 

and Z is equal to zero at that radius which is known as the Schwarzschild radius. Therefore the metric 
(O loses its meaning there. 


4. Internal solutions 

4-1. Distance modulus and redshift 

The distance modulus p and the redshift z of Type la supernovae represent reliable information 
from which we may infer some knowledge about the structure and motion of the universe, see Fig [T| 
That information comes to us at the event (t, r) = (0, 0) through the light of stars which was emitted 
in the past at events (t, rr(t)) with t < 0 and with 0 < rx{t) < R(t). rr{t) is the trajectory of light 
introduced in Section 2.2. /i and z depend on the values of the fields er, to, Z and V at the event 
(t,rr(t)) of emission and for our model that dependence is given by the functions 


z(t) = z(t,r T (t )) = 
Kt) = = 51og 10 



c(t,r T (t)) 


(30) 




These equations are derived in the Appendix, see m and (1551) 9 respectively. Therefore our model 
of the universe if it is valid - must have solutions er (t,r), m(t,r) (or S(t,r)), Z(t,r ), and V(t,r) 
which are consistent with the observed luminosities and redshifts. Such solutions will be found in this 
chapter. 

4-2. Initial conditions a(ti,r) and rh(ti,r) 

We take the initial time as some ti < 0 and our plan is to calculate forward from there into 
the range t > ti and backwards into the range t < ti. Inspection of (l22l) i 9 shows that we need to 
impose initial conditions a(ti,r) and rh(ti,r ) in order to determine Z(ti,r ) and V(ti,r), the former 
by integration of ( 1221 ) 9 . The constant of integration may be found by the assumed continuity of Z at 
the surface R(U) of the initial sphere. But smoothness of Z - also assumed - needs restrictions on 

a ^ ^ ^ Qtti 

the initial functions a(ti,r) and ?h(f,;,r); by ( 1221 b 9 it requires that <r(tj, i?(tj)) and both 
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vanish. Also, obviously, we must have rh(ii, 0) = 0 and rh(ti, R(ti)) = 1 and, for reasons of symmetry, 
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Apart from these conditions a(ti,r) and m(ti,r) are arbitrary or, at least, there is nothing in the 
theory that would restrict their generality. Such a circumstance is commonplace in the solution of 
partial differential equations, the situation merely needs the input of observed - or otherwise given - 
initial data. Here, however, in an application of partial differential equations to the universe, there is 
a problem: Observations are either impossible or insufficient. For instance an observer at the event 
(fi, 0 ) cannot possibly measure a(U,r), because light with its finite speed needs time to travel from r 
to the center at r = 0 . 

In this situation we have to determine the initial conditions a(ti,r) and m(f,,r) in a complex 
argumentative loop as follows. We assume trial functions & tri a i(ti,r), fhtriai (t,; , r) and determine the 
corresponding solution a(t,r), rh(t,r), Z(t,f ), V(t,r ) of the Einstein equations, hence the trajectory 
rxit) and hence the distance modulus /x(f,rV(£)) and the redshift z(t,rTit)) from (1501) and finally 
^l(z) by elimination of t. If the agreement with the observed /r(z)-relation is good, so were the trial 
functions. If not, we have to change the trial functions and try again until we have a satisfactory 
agreement. 

Such a procedure is known as the solution of an inverse problem, a complex variant of the simple 
shooting method. It is a cumbersome and time-consuming process. However, it can be automated to 
a certain extent. Also the process may be abbreviated by an intelligent choice of the trial functions 
and it may be directed by intermediate results. If well-conducted, the process in the end amounts to 
the adjustment of a few parameters to the observed data. 

As specific trial functions for a(ti,r) and m(ti,f) we choose 


f’') — &(tii 0 ) | 1 

dlTltrial 
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dr 
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rn trial (A, f) itself follows by integration of lEfflE over f between 0 and f , since m(tj, 0 ) = 0 holds. 
And a(ti,0) may be calculated from the requirement that m(ti, R(ti)) be equal to 1. Fairly obviously 
the trial functions obey all the conditions we have required. In particular, the form of B(f ) ensures 
that V 2 > 0 holds. This leaves us with five parameters in (l3H) to be determined in the laborious 
adaptive process of solution of our inverse problem, namely n, a, b i, & 3 , 65 . Also Q , the dimensionless 
Schwarzschild radius, and I?(0), the present radius, enter into this problem. Two convenient choices 
of the outer radius R(ii) of the sphere at two initial times ti are R(ti) = 0.44 and R{ti) = 0.25. After 
the calculation the corresponding times ti may be read off from Fig. 3] as ti = —0.677 and ti = —0.789 
respectively, because i?( 0 ) must be equal to 1 . 

We do not exhibit the adaptive process in detail. However, we anticipate its result in order to be 
able to discuss the subsequent steps: As it turns out, the triple Q, n , a need not necessarily be involved 
in the adaptation in order to obtain good results, only the b' s, and i?(0). We find that two choices 
of the parameters, denoted by A and B, give nearly perfect results in the sense that the calculated 
function lies on the observed points in the fashion demonstrated in Fig. [I] These "good" choices 
are listed in Table [T] And the corresponding good initial functions a(ti,f) and are exhibited 

in Fig. [2j 


Table 1: Two "good" parameter choices in the initial conditions for a(ti,r) and at time ti. 



Q 

n 

a 

bi 



R( o) 

A 

0.11 

2 

i 

1.77 

-0.94 

17.65 

10.79 * 10 9 ly 

B 

0.015 

3 

0 

3.9125 

0 

0 

14.27* 10 9 ly 
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Figure 2: Initial functions <r(t;,f) and rh(t;,r) for the parameter choices A (solid) and B (dashed). 


The fact that there are two "good" sets of trial functions emphasizes the non-uniqueness of the 
solution of our inverse problem. There are sure to be many more good solutions, meaning that more 
astronomical observations - in addition to fi- and z-measurements - are needed to describe the universe 
properly. 

4-3. Initial functions for Z(ti, r) and V(U, r) 

Given a (f;, f) and m(ti, r) we may now calculate the initial functions Z(ti 1 r) and V ( ti,r ) from the 
Einstein equations ll221h ■?. Integration of (1^1) 9 provides Z(ti, f) and subsequent insertion of this initial 
function into (j22]h provides V(U,r). The graphs are shown in Fig. [3] 




Figure 3: Initial functions for Z(ti, r) and V(ti,r) for parameter choices A (solid) and B (dashed). The graphs of Z(ti,r) 
combine the internal solution with the external solution (f28l) 


4-4- Step-wise solution for a(t,f), rh(t,r) , V(t,f ), Z(t,r) and R(t) 

We proceed to find the solution of the Einstein equations (122f) . or (E51) . or (1241) . All functions 
depend on (f, f). 

The algebraic equation (l24l) i is only used to calculate the initial condition for V. In order to avoid 
the indeterminate expression (3crf 2 — m')/m' in (l24|h at f = R(t) we use the set (l23l) for the numerical 
solution instead of (1241) . This set represents a coupled system of three partial differential equations for 
m, <r, V and one ordinary differential equation for Z. 

For a very small time step At we replace V in the equation (1231W by the initial condition. The 
result is the simple Burgers equation 


dmif^r) y ~ dmi^r) 
dt + U ’ ’ dr 


= 0 , 


in which m(i,r ) is the only unknown. For the initial condition of m we have rh(ti,r) by (1311) 9 on the 
interval 0 < f < R{tf). The solution follows by the numerical method of characteristics on the interval 
0 < r < R(ti + At). The equation (l23l) i will also give a Burgers equation for a by replacing rh and V 
by their initial conditions. This equation is again solved by the numerical method of characteristics. 
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Figure 4: Rif ) and Vnjt) for parameter choices A (solid) and B (dashed) 




r 


Figure 5: (a): d(t,f). (b): rr(r') along the trajectory. Parameter choice A. 


Now we have rh(t, r) and <f (i, f) and these functions are inserted in the equation (12311 9 for Z , which 
may be solved by a standard method for ordinary differential equations. 

At the end of the first time step we solve the partial differential equation (l23lu for V in the same 
manner as the equation for a. 

The solution rh(i,r), a(i,r), Z(t,f), on tj < t < U + At and 0 < r < R(U + At) is then 

used as initial condition for the next time step. 

4-5. Summary of results 

First and foremost the solution of the Einstein equations provides the radius of the universe as 
a function of time as shown in Fig. 0^,. The present time t = 0 is identified as the abscissa of the 
radius R = 1 . Thus we see that the time U, which we have arbitrarily chosen as the initial time 
with R(ti) = 0.44 for choice A is equal to t, = —0.677, while for R(ti) = 0.25 in choice B it is equal 
U = —0.789. R(t.) grows indefinitely in the future toward an asymptotic slope. Also inspection of 
the graph R(t) in the past reveals an asymptotic approach of R{t) toward the Schwarzschild radius 
Q as t tends to small values. Thus we conclude that the expanding universe has emerged as a sphere 
with that radius. Fig. [4jo shows plots of Vr( t), the surface velocity of the expanding sphere which 
results from differentiation of R(t). It starts with zero velocity at small times and continues to grow 
asymptotically toward Vr = 1. The growth of R{i) and Vr(J:) seems to be accelerating, since the slope 
of V R {t,) is positive and since the graph of R{t) is convex, i.e. it has a positive curvature; at least 
Vr( t) does not stop growing, as far as we can tell, even for large times. We come back to this point in 
Chapter 5. 

The main results of our calculation are the fields a(t,r) , m(t,f) , V^J;) and Z(t,r). Since they 
are fields on the (t, r)-plane, the best way to represent them is by 3-D plotcl- Such plots are shown in 
Fig. [5] through [8] for times —1.5 <t< 0.5 or —2.0 < t < 0.5 and in the interval 0 < r < R(t). The 


8 An alternative is to show movies , rh(t,f ) , Z(i. r) in time t for 0 < r < R(t). That presentation will 

be used later in Chapter 6, at least for cr(t, r) and V(t , r), of which we present screenshots for successive times. 
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Figure 7: (a): Zi n t(t,f). (b): Z(r ) along the trajectory. Parameter choice A. 


plots use different perspectives, depending on whatever seems most appropriate to us. All 3-D plots 
represent solutions for the parameter choice A. 

Additional important information furnished by the Einstein equations are shown in Fig. [31 which 
represent the fields 


S(t, r) and 



— Or 


= C 


\Z{t,r)\ 

S(t,r) 


(32) 


S(t,r) follows from m(t,r) by (TTO1) and c x , the speed of light, is recalled from (O ■ The corresponding 
3-D plots are shown in Fig. [9j 

Fig. [3J:> represents the speed of light at every event (f, r ). Starting from the event (f, r) = (0, 0) we 
may thus determine - by numerical integration - the past trajectory rx{t) of a light ray that reaches 
the center at the present time. That trajectory is shown in Fig. [lOl see also Fig. [22] in the Appendix, 
Chapter 7. It is only from points on this trajectory that we, at (0,0), can receive information about 
light-emitting stars. 

The 3-D plots [5] through [9] exhibit the values of the fields Z(t, rr(t)) through c x (t. rr{t)) along the 
solid lines starting at (t,r) = (0,0)and ending on the surface of the sphere. 

The initial values that start from (t, f) = (—0.677, 0) are also shown as solid lines in the 3-D plots 
5 through 8. Note that those lines in the (<r, t, f)-plot and in the (m,t, r)-plot represent our trial 
functions (15T1) for the “good” parameter choice A. 

The radius R(t) and the trajectory rr(f) intersect for choice A at the event (to, ?o) = (—0.612, 0.489). 
That event is the farthest in time and space which we can observe. It occurred at a time when the 
size of the universe was roughly half its present size; indeed we have tq ps 0.5i?(0). 

Naturally the values of our fields along the trajectory are most important for the astronomer. We 
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Figure 8: (a): V(t,r). (b): V(r) along the trajectory. Parameter choice A. 




Figure 9: (a ):S(t,r). (b): c x (t,r). Parameter choice A 
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Figure 10: ry(£) is the light trajectory. Parameter choice A (solid) and B (dashed). Note that rx(£) is slightly concave, 
see also Fig. 1221 in the Appendix, Chapter 7. 

denote them by a(t) = a(t, frit)) etc., if represented as functions of t with to < t < 0 , or as <j(r) etc., 
if represented as functions r with 0 < r < tq. 

Fig. [5 Jd shows <r(f); it exhibits a non-monotone character due to the competing facts that a was 
larger in the past in an expanding universe but has to drop to zero at the surface of the universe. The 
starting value a ~ 4 means that the central density cr(0, rx{0)) at the present time is four times bigger 
than the mean density at that time, see m*- 

Note that the slope of the light trajectory in Fig. [TO] is the speed of light c x at the event (t, rx(t)). 
It is everywhere slightly smaller than c, but it grows as the central event (0,0) is approached. 

Although the 3-D plots of Figs. [5] through [3] represent the solution of the Einstein equations, it 
is best to postpone a discussion of its qualitative and quantitative features. Indeed, the solution is 
represented here in terms of Schwarzschild coordinates (t,r). These coordinates have allowed us a 
fairly easy solution, but they are bad for interpretation in terms of intuitive notions of space and 
time, particularly time. A case in point is the asymptotic behavior of R(t) at small times, where R 
approaches the Schwarzschild radius asymptotically, whereas we expect an increasingly rapid decrease 
of R toward Q under the action of the gravitatioqj. According to our solution, however, the derivative 
dt/dR is singular for R — > Q and that forces us to deliberate about the rates of clocks in a gravitational 
field as compared to their rates in a Lorentz frame. We postpone this topic because at the present 
stage we do not yet need it; see however Chapter 5. 

What we do discuss next is the question whether our solution, - albeit in Schwarzschild coordinates 
-, does allow us to derive the observed //(z)-relation of Fig. [I]as we have anticipated. Since that relation 
does not explicitly contain r and t, it is valid irrespective of our choice of space-time coordinates. 

4-6. Cosmological conundrum 

In the Appendix, Chapter 7, we have derived expressions for z(t) and n(t), see also (1.'jQ[) . And now 
we have exhibited all functions on the right-hand-sides of these, - namely rr(t), V(t, rr(t)), c x (t, rr(t)), 
S(t,rr(t )) -, so that we may plot z(t) and n(t), or z(r) and n(r). Those plots are shown in Fig. Qj] 
for choice A of the parameters. Elimination of t or r provides the graph fi(z) of Fig. [T] which fits the 
observed data well. This, of course, is a foregone conclusion since choice A represents a “good” choice 
in the sense of Table [I] And the criterion for “goodness” was a good fit of calculated values of the 
n(z)- curve with the measured data. 

Choice B is another good choice and accordingly the n(z) plot for choice B, shown in Fig. [12] is 
again good. In fact, to the naked eye there is no difference between Figs. [T] and [12] 

We conclude that the fi(z) results of Figs. [T] and [12] do not call for a dark-energy-hypothesis nor 
do they require an accelerated expansion as we shall see in Chapter 5. That is satisfactory! But there 
is also non-uniqueness: Although the parameter choices A and B predict the same Hubble diagram, 
they do also predict different density functions a(t ) or a(r) along the trajectories as is illustrated in 
Fig. [13] Other functions, like m[r) or Z(r) will differ as well between the choices A and B. This means 
that we need more observations; observations other than those of /i and z. 


9 Note that the graph R(t) may be read forwards or backwards. In the forward mode, — from left to right — it represents 
expansion and in the backward mode it represents contraction. 


14 









r r 

Figure 11: (a): z(r). (b): and Along the trajectory. Parameter choice A. 
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Figure 12: The fi(z) -plot for the parameter choice B of the universe. 



Figure 13: Density distribution along the trajectory. Choice A (solid), choice B (dashed). 
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For instance, suppose that astronomers had measured crfr) along the trajectory. We should then 
have to perform our inverse adjustment of initial conditions not only for the observed Hubble function 
p(z) as the target function but also for the observed values crfr). However, cr(r)-observations do not 
seem to be available. What comes closest to them are galaxy counts. 

Indeed, among the cosmological observations other than the /i(z)-plot of Fig. |T) galaxy redshift 
surveys [lf|| have furnished a reliable histogram of the fraction of 250.000 galaxies over the redshift 
of the light emitted by them. The histogram has gelled into a simple smooth analytic approximation 
curve of the form (see [lj| p. 1059ff) 

fraction = z 2 exp(—49.3Z 1 ' 55 ). (33) 

This formula may be converted into a "smeared out" galaxy density, and there is the temptation to 
consider this galaxy density as proportional to the matter density a. However, we are not certain 
whether far-away galaxies have the same mass as close ones. Therefore we have not incorporated 
galaxy counts - reliable as they may be - into our adaptive inverse scheme. That remains to be done 
in the future after more deliberation. 

Moreover, it can hardly be expected that an eventual observation of crfr) would suffice to identify 
unique initial conditions. More observational evidence is likely to be needed and the question is: Which 
additional observations are feasible and what solutions do they entail? That question represents the 
cosmological conundrum. 


5. Interpretation of results 

5.1. Radius of the universe as a function of the proper time r=0(t,0) at the center r=0. 

As was mentioned before, - toward the end of Section 4.5 -, the use of Schwarzschild coordinates 
for the solution of the Einstein equations has made it difficult, or impossible, to interpret the results 
in terms of our intuitive notion of time. We believe that it is the proper time 

r = Oft, 0) 


of the observer at the center of the universe that comes closest to our intuitive notion of time and that 


Vt (t,r) 


dr\ Vft,r) 

^) e = w 


(34) 


corresponds to our intuitive notion of velocity. Therefore we proceed to replace t by rft). 

d0 

From © and the knowledge of V(t,r), Z(t,r), and S(t,r), see Figs. 1817191 we obtain —( t,r ) and, 
d dO dt ' 

in particular, the graph — ft) = —(t, 0) of Fig. [TT] Integration over t gives rft) to within a constant 
dt. dt 

which we identify by requiring that r = 0 holds for t = 0. The resulting graph rft) is also shown in 
Fig. [TT] Inspection shows that for large values of t the proper time r becomes asymptotically equal 
to t , while for small values of t the proper time becomes independent of t. 

We conclude that sufficiently far in the past the tiniest increase of r requires a large increase in t, 
or else: the rate of Schwarzschild coordinate clocks, which measure t, is much larger than the rate of a 
clock at rest in a Lorentz frame at the center. It is for this reason that the radius R(t) in Figs. [TUI and 
M does not seem to change at all with time t as the radius approaches the Schwarzschild radius Q. 

R(t) on the other hand - obtained by elimination of t between R(t) and r = Oft, 0) - seems to 
grow linearly starting with R = Q, see Fig. [T5k which shows both functions Rft) and R(t). Closer 
inspection, however, reveals that R(r) is a slightly concave function so that the observer at r = 0 
sees a decelerating expansion in terms of his time r. Fig. [T~5b demonstrates this fact by focussing the 
attention on the neighborhood of the point where R(t) emerges from the Schwarzschild radius; that 
happens at approximately r = —0.85. The fact is further illustrated in Fig. [TU] where the velocity of 
expansion dRfr)/dr is plotted. The figure shows that the velocity is a convex function and has not 
even developed into a straight line at r = 0, i.e. at the present time. For large values of r the surface 
velocity tends to a constant value < 1. For choice B this is obvious from inspection of Fig. [TUI for 
choice A we need to pursue the graph far to right in order to see it drop below the value 1. 

The plots of Fig. [TT] correspond to those of Fig. IT5l except that they are calculated for the parameter 
choice B. For both choices, A and B, we see the same qualitative behaviour of the i?(r)-curve. The 
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Figure 17: Parameter choice B. Radii R(t) and R(r). Also trajectories rj-(t) and rj-(r). 


concavity in that curve - as exhibited in Fig. 117b - is a little more pronounced for choice B than for 
choice A. From both figures we conclude that the expansion of the universe is decelerating, since R(f) 
is concave. 

The significant difference between R(t) and R(r) in Figs. 1151 and 1171 - one convex and the other concave - is 
reminiscent of the case of a mass which drops into a black hole. That case is routinely treated in books on general 
relativity, e.g. see Shapiro and Teukolsky There too the distance h(t) of the mass from the black hole tends to a 

finite value as the Schwarzschild radius of the mass is approached, while h(r), where r is the proper time of the mass, 
drops precipitously toward zero. 


Our solution does not provide graphs for R(t) < Q simply because we do not get values for 
R(t) < Q. What we do obtain is the proper time when the expanding sphere emerges from the 
Schwarzschild radius. That time may be read off from the figures as the abscissae of the end points of 
the R(t)- curve and maybe that time is appropriately called the "time of emergence." We know nothing 
about the period of time before that age, which was where the "big bang" happened - if it happened 
- and where the expansion began. Our cosmological model is not able to provide information about 
such phenomena, at least not in the present form. 

5.2. Trajectories 

The Figs. US and M also show graphs of the trajectories rxit) and of tt{t ); the latter results by 
elimination of t between rx(t) and r(t). It is only from events on these curves that we may receive 
information by light about the past. In analogy to Fig. [TO] we denote the coordinates of the points 
of intersection of the graphs tt(t) and R(t) by (to,to). We cannot observe what happened before to 
and at a greater distance than rg. Therefore to may be called a " look back time." Its value is given in 
Table [2] along with other interesting properties of our models. 

While the trajectory rr{t) is very slightly concave, the graph Vt{t) is straight and its slope is equal 
to 1 which, in dimensional terms, means that - of course - the light moving along tt (t) has the speed 
c of light in the Lorentz frame at the center. 

5.3. Some results in tabular form 

We summarize some of our results in Table [5] We recall that I?(0) was also listed in Table |T] 
and that it was determined by adjusting our parameters to give a good agreement between theory 
and observation of the /i(z)-curve. Given 7?.(0) the total mass follows from the definition (EB1) of the 
Schwarzschild radius. 


Table 2: Some results for the parameter choice A and B 



R( o) 

M 

Time of emergence 

TO 

ro 

Ho 

A 

10.79 Gly 

1.5 14- 10 52 kg 

« -0.85 = —9.22Gy 

= —5.311Gy 

= 5.281Gly 

0.22 • 10 _17 s -1 

B 

14.27 Gly 

0.273 • 10 52 kg 

« -1.00 “ -14.26Gy 

“ —7.212Gy 

“ 7.210Gly 

0.22- lO^V 1 
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Figure 18: Velocity as function of distance along the trajectory. Choice A (solid), B (dashed). The slope of the tangent 
may be identified with the Hubble constant. 


The data given in the table for choices A and B differ widely. And, what is more, they differ 
considerably from the values provided by the currently popular FRW cosmology, where, for instance, 
the age of the universe is confidently given as —13.7Gy. Well, we can only say that that value is heavily 
dependent on the Robertson-Walker model along with the concept of a hypothetical dark energy which, 
moreover, makes up 70% (!) of all mass in the universe. We do not consider that notion convincing. 


5-4■ Hubble's law 

The very fact that the function t = t{r ), according to Fig. [ll] has a singularity for r « —0.85 
disqualifies t as an appropriate measure for time. That fact makes it impossible to interpret V(t,r ) = 

(t, r ) as a reasonable measurable velocity. Therefore we ought to reformulate the calculated 

e . 

field V(t,r) in terms of the proper time of the observer in the center. As in (1M1) . we consider 



V T (t,r) = 


V(t,r) 

w(t) 


as the proper measure of velocity. This is a function on the (t, r)-plane. From it the velocity V T (t, rj'(i)) 
along the trajectory, i.e. the velocity of a star which has sent out light at a time t < 0, may be 
calculated. Hence follows the velocity V T (r) of the star at the distance r along the trajectory by 
elimination of t between V T (t,rT(t )) and ry(f) ■ It is shown in Fig. [151 and we see that it increases 
with an increasing gradient as r approaches the surface of the sphere so that on the surface - where 
r = ro - the velocity V T (r) is close to 1. This means that the surface moves away from us roughly with 
the speed of light. 

Close to the center, where Hubble Q made his observations, the convex curve V T (r) may be 
approximated by its tangent whose slope is the Hubble constant Hq = 0.759 for choice A, see Fig 


1181 In dimensional form this value is thus given by Hq = Hq 


R( 0) 


= 0.22 


10 _17 s _1 ; 


it is listed 


Table [2] This value agrees well with the value Hq = 0.226 • 10 _17 s _1 which is most often given in the 
astronomical literature, e.g. see [lj. 

For choice B the initial slope of the V T (r) curve is given by Hq = 0.986 which again implies 
Hq = 0.22 • 10~ 17 s _1 in dimensional form, see Table [2] We must not be surprised that the choices A 
and B agree as to the value of the Hubble constant despite their disagreement about all other values 
in Table [2] Because indeed, the models do agree on the p(z)-relation, see Figs. [Tl and l~H?l and Hubble's 
law is a corollary of that relation. We proceed to describe this point briefly. 


Generally we have, by 455)1 in the Appendix, Chapter 7 

m(z) 

dL = 10 pc 10 5 , 

where d l is the luminosity distance defined in 454). For small values of z we have from the observed data, see Figs. []] 
and 1121 

fi(z) rj 33.38 + 5 log 10 -. 

v ' 10 0.01 
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Figure 19: (a) Newtonian deceleration, (b) Repulsive acceleration, (c) Terms with Z. (d) Overall deceleration along the 
trajectory. 


Also, by {54} i and {S3 we have for small z 

cLl r and z V T . 

Hence follows for the non-dimensional Hubble constant 

- Vr R( 0) 1 

0 f 10pc 10 WP' 

c ~ 

For the dimensional Hubble constant H o = ^ Ho the choice-dependent factor -R(O) cancels and thus Ho is the same 
for choices A and B: Hubble’s observations for small z have no independent status vis-a-vis the /r(z) observations. 


5.5. Accelerating and decelerating contributions to the expansion of the universe 


We have argued that the Schwarzschild coordinate time t is not a good measure of time and that 
the proper time r = 9{t 1 0) in the center of the sphere is better. The difference is not trivial; indeed 
R(t) exhibits an accelerated growth while R{r) exhibits a decelerated growth, see Section 5.1. To be 
sure, according to Figs. [15] and [l7] both modes of growth eventually settle into a uniform expansion 
with a constant velocity close to that of light, see Fig. [TH] This is strange in itself, quite apart from 
the question whether the radius in the past has been accelerated or decelerated toward that uniform 
value. Indeed, do we not expect - as more or less classical physicists - that there should be a continued 
deceleration under the effect of the gravitational attraction? In answer to this question let us discuss 
the role of gravity in our model; that role has so far been disguised by the unfamiliar form of the 
equations and by the fog raised through the numerical solution of our inverse problem. 

We call the attention to (l22lh which, after a little rearrangement and the insertion of dimensionless 
quantities reads 


dV_ 

dt 


-\q% 


Q^v 2 


1 -Q 1 


1 to IV dZ 

-Q—(Z - 1) -t-- 

2 r 2V ; 2 Z dt 


(35) 


in terms of the Schwarzschild time an 




dVr _ l p m | pda- 2 1 1 to / Z(r,r) _ \ 1 - / 1 dZ{r, r) _ 1 dZ(r,0) \ 

dr 2 y r 2+y r 2 T 1 - Qf 2^r 2 [z(r,0) J 2 T [z(r,r) dr Z(r, 0) dr ) 

(36) 

in terms of the proper time r(t). 

It is clear that those equations must be identities when we insert the solution exhibited in Chapter 4, 
because the solution has been derived from them. However, (1351) and (1361) may be viewed as expressions 
for the accelerations dV/dt or dV T /dr. We concentrate on dV T /dr , because that is the acceleration as 
seen by the observer in the center. The first term on the right-hand-side is the classical gravitational 
acceleration. It is clearly attractive because of the minus sign. But it is not alone! 


10 The hats for non-dimensional quantities are now dropped. Besides the tilde denote functions of r as the time variable. 
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Fig. [Tin shows plots of the four contributions on the right hand side of (1551) along the trajectory 
r-r(t): In (a) we see the attractive classical (Newtonian) contribution, while (b) shows the repulsive 
second term which grows with growing velocity as the surface is approached. In (c) we have plotted 
the third and fourth terms, - the terms with Z and its derivatives; those are negative and therefore 
attractive. Finally the fat graph (d) exhibits the entire right-hand-side of (1361) . We see that the 
classical attraction is diminished by the other terms throughout the whole spread of the universe but 
that it remains attractive, albeit weakly near the surface. We conclude that there is no accelerated 
expansion anywhere. Therefore there is no need to introduce a hypothetical dark energy to create the 
accelerated expansion. 

Also we now understand why it is that the surface of the universe moves with a nearly constant 
speed for large times despite of what we are accustomed to consider as the gravitational pull. Indeed, 
the gravitational pull is quite small on the surface. 

The parameters for which the graphs of Fig. Q12] are drawn, are those of parameter choice A. For 
choice B we obtain smaller values - obviously reflecting the smaller mass M of choice B - and the final 
value at f = —0.5 remains negative: Apparently the deceleration has not come to an end yet at the 
surface of the sphere. 

The various accelerations of Fig. [Til] are all calculated for the present trajectory, i.e. the trajectory 
that passes through the event (0,0). Surely the plots will look different for trajectories at earlier times. 
And it is even conceivable that for such earlier times the graph (d) may exhibit an overall acceleration 
instead of an overall deceleration. In that respect it seems worthy of note that for real early trajectories 
Z = 0 holds and that the accelerating term (b) in Fig. Qj|]becomes singular, see Fig. [TJi and eqn. (1551) . 
That aspect, which touches on the notion of an inflationary past period of the universe, will be the 
subject of a subsequent study, if the present one fares well. 

In the present study the main conclusion is that there is no overall acceleration along the present 
trajectory and therefore there is no need for a hypothetical dark energy. 

6. Remark on the cosmological principle 

We are fully aware of the fact that our model of a sherical universe floating in infinite empty space 
violates the cosmological principle according to which we, the observers should not occupy a privileged 
place in the cosmos. In our model clearly the center of the sphere is privileged and yet that site is 
assumed to be the place of the observer. 

In early cosmology the problem became urgent with Hubble’s observation that the galaxies in our 
neighbourhood move away from us with a speed proportional to their distance. The question was: 
Why are we thus privileged? That question lost some of its urgency when it was recognized that in an 
infinitely extended homogeneously dense matter distribution with an expansive velocity proportional 
to the distance from one point the matter is in fact moving away from all points in the same manner. 
So, in such a universe Hubble’s observation does not put us in a privileged spot. 

Of course, if the universe is spherical with a finite radius and not homogeneously filled with matter 
- as in our case - the argument does not hold, or at least it does not hold exactly and everywhere. 
We may conjecture, however, that the argument does hold approximately in an inner sphere around 
the center where the density is nearly constant and where the matter moves away from the center in 
the Hubble way. Fig. l20l and [2T1 show that there is such an inner sphere, well pronounced for choice 
B and somewhat less well pronounced for choice A. 

Clearly that argument needs strengthening. It is offered here loosely so as to anticipate the objection 
that our model of the universe places the observer in the central point of the universe. We believe that 
the model is still good, if the observer is placed anywhere within the inner sphere of near-homogeneous 
density. 

7. Appendix: Redshift of Type la supernovae and luminosities. 

7.1. Scope 

The study of redshifts requires an investigation of the Doppler shift and aberration and of the 
gravitational frequency shift. We discuss these phenomena in the subsequent Sections 7.2 through 7.4 
with (1511) and Fig. fTTh as the result. Section 7.5 is given to a discussion of the apparent luminosity 
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Figure 20: a(t, r) for 0 < r < R(t) at various times, (a): Choice A. (b): Choice B. 




Figure 21: V(t,r) for 0 < r < R(t) at various times, (a): Choice A. (b): Choice B. 


with (|55lt > and Fig. [Tib as the result. Both results (IFTl) and (|5R|) o have been anticipated in Sections 
4.1 and 4.6. 

The formulae for redshifts and luminosities differ between different cosmological models, because 
the concepts of frequency and wave length are non-trivial in relativity. Etherington [l9| has studied the 
problem as early as 1933. See also Ellis [2C| who reviews Etherington's work as a "Golden Oldie." Still, 
however, there are conflicting formulae used as relations between absolute and apparent luminosities 
in the modern literature. Thus the relevant results of the Friedmann-Robertson-Walker model are 
useless for our more mundane model; there is not even a clear-cut Doppler effect in the FRW-model, 
see [lj. Therefore we have rederived the relevant formula in this chapter and we hope and trust that 
we got things right. 

7.2. Dopplershift. and aberration 

In this chapter we consider light emitted by a star or galaxy moving with the radial velocity V 
with respect to the center of the sphere and absorbed in the center. 

We look at the light from three different frames of reference: i.) the local and momentary Lorentz 
frames z a and z^ - with coordinates 9 , z 1 , z 2 , z 3 and 9q , Zq, Zq, Zq respectively - which accompany the 
star and the observer and ii.) the frame x a - with coordinates t , ar, x 2 , a: 3 - in which we have solved 
the Einstein equations. It is now appropriate to introduce rectangular Cartesian spatial coordinates 
z a and x“, since the light does not necessarily move in the radial direction. However, we let z 1 and 
a: 1 point into the radial direction so that the unit vectors of propagation of the light may be given in 
terms of the direction angles a and f3 by 

n a = ( cos a, sin a cos /3, sin a sin ft I and n a = ( cos a, sin a cos /3, sin a sin/? I . (37) 

z \ z z z z z) x \ x x x X z) 

Let the light locally be described as a plane wave with frequency v and wave length A. The space-time 
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wave vector k a is then given by 


z 1 \ 

k a = 27t —, ——n a , where vA = c 

z \ c Az I zz 


k a = 2-k\ —, ——n a , where = c. 
x \ c Ax ] 


XX X 


( 38 ) 


We consider k a as given and calculate k a using the transformation matrix between the frames {z a } 
and {/} 

(39) 


k ?~ dx« k t 

The transformation matrix may be calculated from the invariance of the infinitesimal distance element 
in space-time, see ®, © 

ds 2 = — c 2 dd 2 + dz 12 + dz 22 + dz 32 = Zc 2 dt 2 + Sdx 12 + Sdx 22 + Sdx 32 . (40) 

The calculation provides 



( V\Z\ 

V V\Z\ 

0 

0 \ 




dz p 

v Vs 

Vs 

0 

0 

dx a ~ 


fv 


0 

0 


0 


V o 

0 

0 

/ 


(41) 


f dv 

As before, V = — ) is the velocity of the star which emits the light; it is non-negative for the 


\dt 

expanding sphere. Insertion of (l4lT) into (l39l) gives 
v 1 + T-ri 1 




=, ~ nl = 'M 


n 1 + ~ v 


1 - 


v 2 V c 


1 - 


v 2 V 


, -n 2 = V\Z\n 2 , - n 3 = \/\Z\n 3 

z y x 


(42) 


Hence follows in terms of the direction angles a and j3 of the light 


- = V\z\- 


1 + — cos a 


cos a = 


1 - 


y 2 


cos a + 

Z X 

1 + — cos a 


,P = (3. 

X Z 


(43) 


These are the equations of Doppler shift and aberration, so called in analogy to similar effects in 
acoustics. The aberration equations (1431) 9 ^ imply for an element dkl = sin adad/3 covered by light rays 
of neighbouring directions 


i v 
cos a + — 


dfl = -tt- 

* 1 + — cos a 


1 - 


V 


1 + — cos a c 

~ Z X 


dfi. 


(44) 


In particular, when the light moves radially inwards, i.e. for a = 7r, we obtain 

v 

£ = y/\z\ - -and dn = - - ^dn 


\ 


1 + f 


l- A 


(45) 


so that the Doppler shift is a red shift in this case and the solid angle element is bigger in frame {x a } 
than in the Lorentz frame {z a }. 
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Figure 22: Two trajectories of light emitted at (t, r) (solid) and (t + A t e , r) (dashed) respectively and absorbed at (0, 0) 
and (At a , 0). The black bars demonstrate that the horizontal distance of the graphs has become smaller when the light 
reaches the origin at r = 0.Parameter choice A. 


7. 3. Gravitational shift 


The Doppler redshift is not the only phenomenon that affects the frequency of the light emitted 
by a star at (t,r) and received by an observer at (0,0). There is also a blueshift due to gravitation, 
because, after all, the light gains energy by "falling" toward the center. In order to describe that 
additional frequency shift, let us look at the trajectory r'r(f) again, the curve exhibited in previous 
figures and reproduced again in Fig. [22] The light - already redshifted from v to v according to flSlb 


- starts at (t, r) and proceeds along the lower graph in Fig. [22]to the origin at r = 0. Now let there be 
a second trajectory of a light ray emitted at the same point but at the later time t + A t e and absorbed 
at the origin at time A t a . 

Obviously, since the trajectories are concave and roughly "parallel", we have A t a < A t e . The short 
solid bars in Fig. [22] illustrate and emphasize the situation. A moment's reflection shows that the 
A t e 

ratio tends to the inverse of the slopes of the trajectories as At e tends to zero: 


A t e 

A t a At e —>0 


d,TT I 

dt IQ 

drT I 
dt 1 1 


c( 0,0) 

X 

c(t,r T (t))' 

X 


(46) 


Now, let the two emissions be consecutive emissions of the maxima of a harmonic light wave so that 

A t e = — and A t a = — hold, where v is the frequency of the light at (0, 0). In that case we have 
v v xo 

x xQ 


V 
xO 

V 

x 


c( 0,0) 

X 

c(t,r T {t))' 

X 


(47) 


This means that the light arriving at the origin has been blueshifted, since the right hand side of 
m i s bigger than 1, see Fig. [22| 

In a final step we ask for the frequency v of the light in the Lorentz frame at the event (0,0) where 

_ zO 

the velocity is zero. In analogy to (05]) i we have 


v 

xO 

V 
z 0 


V\WM, 


x 

xO 


zO 


1 

n/sM' 


(48) 
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Figure 23: The contributions of Doppler shift and gravitational shift and the overall red shift. Note that the gravitational 
contribution is small compared to the Doppler contribution; the graphs refer to choice A. 


7.^. Summary on redshift 

Elimination of v and v between (1E51) . (1471) . (148f) i provides the overall redshift formula 

xO x 


V 
zO 

V 

z 


c(0,0) 

X 

c{t,r T { 0)) 

X 

1 

! + 1 2 G 


\ Z (t, r T(t))\ 

\Z(0, 0)1 ^ 


1 - 


1 + 


1 

i + ZD 


V{t, r T(t)) 

c(t,r T (t)) 

V(t,r T (t.)) 

c(t,r T {t)) 


(49) 


It is customary to introduce a shift factor 2 to replace the frequency quotients such that z is positive for 
a redshift and negative for a blueshift. The definition of 2 is indicated in (H^l) both for the gravitational 
shift and for the Doppler shift. The overall shift factor 2 follows from 


1 + 2 1 + Zq 1 + Zd 

Hence follows z(t, rrit)) for stars on the trajectory as a function of t. Fig. IE51 shows graphs of the two 
contributions to frequency shift and the overall effect which for our choice of parameters is a redshift. 
Combining (l49l) . (1501) and the definition of c x from © we obtain 


2(f) 


I 5(0,0) 
S(t, r T (t)) 


N 


1 1 V(t,r T (t)) 
c(t,r T (t)) 

1 _ V(t,r T (t)) 
c(t,r T (t)) 


- 1. 


(51) 


It follows that the measurement of the shift factor 2 is not equivalent to the measurement of the velocity 
of the emitting star. Indeed, since c x depends on the metric components S(t,rT(t)) and 
those fields influence the relatioiFL 


7.5. Type la supernovae, their luminosities and distances. 

Astrophysicists believe that Type la supernovae represent the collapse of white dwarfs when they 
have accumulated more mass from their neighbourhood than they can carry according to the Chan¬ 
drasekhar limit. It is then plausible to assume that these supernovae all have the same absolute 


1 For V = 0 equation JED and m combine to give the purely gravitational blueshift 2 : = \ 1 — Q —— 1. And without 


11 + V 

gravitation the equations reduce to the relativistic Doppler shift z = \/ - --r- — 1. 
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luminosity L , i.e. rate of energy emission. And from observations of supernovae of known distance 
that absolute luminosity has the value L = 7.5 * 10 20 J/s corresponding to an absolute magnitude 
of -190 Let the emission occur at time t at the distance rx(t) on the trajectory of the light. We 
proceed to calculate the apparent luminosity l, i.e. the transmission rate of energy per unit area at 
(t,r) = (0,0), i.e. in the center of the sphere which represents the universe in our model. 

Let dN be a number of photons of energy hv emitted by the star at the distance rr(t) from the 
center into the element of solid angle dfh Because of the isotropy of the emission in the Lorentz frame 

Z 

dn 

z a , we have dN = - -—, where E is the total energy emission. The same number of photons - now 

hv 47T 

er 2 dtt 

with energy hv must pass through the area r 2 dQ in the center of the sphere so that dN = —— — 

ZO x hv 

zo 

holds, where e is the energy transmission per unit area. Hence follows 


^ dfl er 2 dfl 

J^__z_ _ _X 

hv 4-7T hv 

z ZO 


(52) 


The corresponding rates of emission and transmission are 


L = 


—— and l 

At z 


e 


At z 


Elimination of E and e gives 


l = 


L At, 


v dil 


Anr 2 Atzn v dfl 


or with 


At z 

A U, 


v 

ZO 

V 

z 


l = 


L 

Anr 2 



(53) 


This is the desired relation between the apparent luminosity l which is measurable and the absolute 
luminosity L which is known, see above. 

We reformulate the right hand side of (l53l> in terms of the redshift factor z: 

By (H5H and (HiH) through (1511) we have 


or, in abbreviated form 


L 5(0,0) 1 

47rr T (t) 2 S(t, r T (t)) (1 + z) 4 ’ 


l = 


L 

4 


where d l 


r r{t) ] 


I S(t,r T (t)) 
5(0,0) 


(l + ~~) 2 


(54) 


di, defined in (1541) 9 will be called the luminosity distance appropriate for our model. 

In terms of apparent mag nitude m and absolute magnitude M - the luminosity measures preferred 
by astronomers - we hava 13 l 


H = m — M = 5 log 


10 


d_L 

lOpc 


or = 51og 10 


lOpc 


M*)i 


I S(t,r T (t)) 
5(0,0) 


(l + 2 ) 2 


(55) 


12 For a qualification of that categoric statement and secondary effects that may modify it we refer the reader to the 
article fl3| by S. Perlmutter. 

13 We have to accommodate astronomers, because they report their observations - like those of Fig. [T] in a (/r, z)- 
diagram and they use parsec (pc) as a standard distance. We trust that the reader will not confuse the present m and 
M with the partial mass and the total mass which are denoted by m and M in the paper elsewhere. 
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